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We discuss the effect of Kaluza-Klein (KK) modes of bulk metric perturbations on the second 
Randall- Sundrum (RS II) type brane cosmology, taking the possible backreaction in the bulk and 
on the brane into account. KK gravitons may be produced via quantum fluctuations during a de 
Sitter (dS) inflating phase of our brane universe. In an effective 4-dimensional theory in which 
one integrates out the extra-dimensional dependence in the action, KK gravitons are equivalent to 
massive gravitons on the brane with masses m > 3H/2, where H represents the expansion rate of a 
dS brane. Thus production of even a tiny amount of KK gravitons may eventually have a significant 
impact on the late-time brane cosmology. As a first step to quantify the effect of KK gravitons on 
the brane, we calculate the effective energy density and pressure for a single KK mode. Surprisingly, 
we find that a KK mode behaves as cosmic dust with a negative energy density on the brane. We 
note that the bulk energy density of a KK mode is positive definite and there occurs no singular 
phenomenon in the bulk. 

PACS numbers: 04.50.-|-h; 98.80.Cq 



I. INTRODUCTION 



The idea that our Universe might be a brane embedded in a higher dimensional bulk spacetime has attracted 
tremendous attention in the last few years. A particularly attractive framework, especially from a gravitational point 
of view, is the so-called second Randall-Sundrum (RS II) scenario, where our brane- universe is embedded in an Anti-de 
Sitter (AdS) five-dimensional bulk spacetime 1]. If the brane is endowed with a (positive) tension, tuned with respect 
to the (negative) bulk cosmological constant A5 = —6/1'^, then, as shown by Randall and Sundrum, the geometry on 
the brane is Minkowski and the gravity felt on the brane is similar to standard 4D gravity on scales r ^ i Q, |^ . 

A considerable amount of work has also been devoted to the cosmological extension of the RSII model, in order 
to describe the cosmological behavior of a brane in this framework (see 0, 0, IS reviews) . A significant result 
in this perspective has been the realization that the Friedmann equations must be modified when the cosmological 
energy density becomes of the order of the brane tension. Another difference with the standard Friedmann equation is 
the presence of an additional term, usually called dark radiation, which represents the influence of the bulk geometry 
on the brane cosmological evolution. The dark radiation term can be related, from the bulk point of view, to a 
five-dimensional gravitational mass. 

All the fundamental results just discussed have been obtained by assuming from the start an exact cosmological 
symmetry, by which we mean that the bulk spacetime is supposed to be foliated by homogeneous and isotropic three- 
surfaces. At all instants in its history the brane is supposed to coincide with one of these symmetric three-surfaces. 
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However, things are more complicated if one considers a more realistic framework where perturbative deviations from 
homogeneity and isotropy are allowed. Such perturbations must be taken into account, for example, if one wishes to 
confront the predictions of brane cosmology with the high-precision measurements of the CMB anisotropies. 

Even if one is interested only in homogeneous and isotropic cosmology, the existence of perturbations can affect the 
homogeneous evolution. An interesting example is the question of dark radiation. The presence of matter fluctuations 
on the brane generates perturbations of the spacetime metric in the bulk. This process may be regarded as the 
emission of bulk gravitons from the brane. These bulk gravitons will on average contribute to the gravitational mass 
in the bulk, and hence modify the evolution of the dark radiation on the brane. That is, it no longer behaves like a 
free, conserved radiation as in the strictly symmetric case. Detailed analytical and numerical calculations of the effect 
of the bulk graviton generation were performed for a radiation-dominated brane universe 0, 0, ITol ITTl IT^ IT^ [T3 | . 

In the present work, we consider a similar, but slightly different problem. We study the impact of the bulk metric 
perturbations which are generated in the bulk (or present from the beginning) on the homogeneous evolution of the 
brane. The bulk metric perturbations are naturally produced via quantum fluctuations during brane inflation |15| . 
and we will concentrate in this work mainly on a de Sitter (dS) brane. From the 4-dimensional point of view, the 
bulk metric perturbations can be decomposed into a massless zero mode and an infinite number of Kaluza-Klein (KK) 
modes with effective mass m > 3H/2, where H is the expansion rate of the dS brane. The time evolution of the zero 
mode is similar to the standard four dimensional perturbations although its amplitude, as determined by the vacuum 
quantum fluctuations, depends on the energy scale of the dS brane expansion. In contrast, the squared-amplitude 
of KK modes on the brane decays as a^^ and thus becomes rapidly negligible during brane inflation. However, 
after brane inflation the background energy density in a radiation-dominated Friedmann-Lemaitre-Robertson- Walker 
(FLRW) era decays as a~*, hence the massive modes of gravitons may affect the late-time cosmological evolution 
of the brane. Here it should be noted that the concept of KK modes, which assumes the separation of variables 
with respect to the fifth dimensional coordinate, is only approximately defined in general. It is an important but 
longstanding problem to quantify the effect of these approximately defined KK modes on the brane evolution. 

In this paper, in order to discuss the cosmological impact of KK gravitons on the brane cosmology, we derive the 
effective stress-energy of a KK mode on a separable (e.g., a dS brane) background, taking possible backreaction in the 
bulk and on the brane into account. Then we extrapolate our result to a FLRW cosmological background on which 
a KK mode can be approximately defined. We show that a sufficiently massive KK mode, which may constitute a 
non-negligible, if not dominant, fraction of the contribution of all the KK mode, when they are summed up, behaves 
as cosmic dust, which is consistent with the linear perturbations, but the effective energy density is negative. 

This paper is organized as follows. In Section II, we discuss the case of a massless, minimally coupled scalar field 
because the situation is similar to the tensor case but simpler. We find that a massive KK mode behaves like dust 
with negative energy density. In Section HI, we turn to the main topic of this paper, namely, the backreaction of the 
KK gravitons on the cosmology of the brane. We find again that a KK graviton mode behaves as negative energy 
dust. In Section IV, we consider our results from the bulk point of view, and discuss its impact on the cosmological 
evolution. In Section V, we summarize our results. Some useful formulae are given in two appendices. In Appendix A, 
the components of the bulk curvature tensor up to second order in the metric perturbations are given. In Appendix B, 
the computational rules for averaging tensor components that are quadratic in the metric perturbation are given. 



II. THE CASE OF A BULK SCALAR FIELD 



Before tackling the main subject of this paper, that is the backreaction of KK gravitons, it is instructive to discuss 
the case of a massless, minimally coupled scalar field, because the behavior of its perturbations is quite similar to the 
KK gravitons but it is much simpler to analyze [l6| . 

We thus consider a homogeneous scalar field and assume its amplitude (j) to be small so that its effect can be treated 
perturbatively: in particular, the backreaction of the scalar field on the metric will be of order 0{(jP). 

Eventually, we would like to discuss the backreaction for a general cosmological background. However, for the 
general case, it turns out that the field equation (either for the scalar field or, later, for the gravitons) is not separable 
and the notion of a KK mode cannot be well defined. The separability property is satisfied only for two limiting cases. 
One is the case of a de Sitter brane. In this case, the brane is exponentially expanding with a constant Hubble rate 
H and one finds a mass gap Am = SH/2 between the zero mode and KK modes. Thus the continuum of KK modes 
starts above the mass 3H/2. The other case is a low energy cosmological brane, in which case the dependence on the 
extra dimension can be approximated by the profile obtained for a static brane, i.e. the RS brane. 
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A. Einstein scalar theory in the bulk 



We start from the five-dimensional action which consists of the Einstein-Hilbert term, a cosmological constant A5 
and a bulk scalar field, complemented by the four-dimensional action for the brane: 

S=J^f d^xy^( - 2A 



24 



- 2A5) + J d''x^(^-^g''''dacl)db<l> - V{^)) + J d^xV^(^~a + , (2.1) 



where q is the determinant of the induced metric on the brane, which we denote by qaf3, and Cm is the Lagrangian 
density of the matter confined on the brane. The Latin indices {a, &,•■■} and the Greek indices {a, /3, • • •} are used 
for tensors defined in the bulk and on the brane, respectively. We will assume that the brane tension on the brane is 
tuned to its RS value so that K|cr^ = — 6A5. We also take a constant bulk potential 

!/((/.) = > , (2.2) 

so that the scalar field is effectively massless. 

We consider backgrounds given by a fixed value of the scalar field which we choose ((> = 0. For a non-zero Vq, one 
has a de Sitter brane background, which will be discussed in subsection B below. For Vb = 0, one has a low energy 
cosmological brane, discussed in subsection C. 

The field equation for the bulk scalar field is linear and given by 

□ 5(/' = 0. (2.3) 

Since we consider a background configuration with (jj = 0, the solution of the above equation can be seen as a 
perturbation. This perturbation will induce a bulk energy-momentum tensor, of order 0{(j>^), which embodies the 
backreaction of the scalar field on the metric. This is the effect we wish to calculate explicitly. 
The variation of the action H2.1|l yields the five-dimensional Einstein equations 

(^^Gab + A^gab = -K^^agab + ^I'^ab + {-O-Qab + Tab)S{y - J/o) (2.4) 

where we have implicitly assumed a coordinate system in which the brane stays at a fixed location y = yo and where 
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^{V—qCm) (2.5) 



represents the energy-momentum tensor of matter confined on the brane. The stress energy tensor of the bulk scalar 
field, not including the constant potential Vb, is given by 

Tab = 4>,a4',b - -^gabg'"''4>,c4'.d ■ (2.6) 

It is useful to consider the projection of the gravitational equations on the brane^^- Taking into account the bulk 
energy-momentum tensor, one finds 

^'^Gy = -l^VoS'^^ + l^'^r^p + kItC'^^ - , (2.7) 



where 



<^'"<^,/5 + S'^p (1^% - ^<"0,p0,.)] , (2.8) 



and Eafs is the projection on the brane of the bulk Weyl tensor and is traceless by construction. If, in addition, 
one assumes the brane geometry to be homogeneous and isotropic then the components of Eap ( in an appropriate 
coordinate system) reduce to E^t and 



E'j = -ls',E\. (2.9) 

By using the four-dimensional Bianchi identities, and assuming that the brane matter content is conserved, one is 
able to express the component E*t in terms of the values on the brane of the bulk scalar field and its derivatives [Tif : 

E\ = 4 /* dt'a^ ( 9tT(^)S + 3^T(^)S - ^t(''^\) , (2.10) 
7to V a a J 
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B. KK mode on a de Sitter brane 



First, we consider the case of a de Sitter brane. The bulk metric around a de Sitter brane can be expressed as 

ds^ = dy^ + l?{y)-^^^dx>'dx'' , (2.11) 

where the warp factor b{y) is given by 

biy)=H£smHy/£), (2.12) 
and 7p,y is the 4-dimensional de Sitter metric, which may be expressed by using a flat shcing for simphcity: 

^^vdx^^dx" = -dt^ + a'^{t)Sijdx'dx^ , 

a{t) = e"' , = ^kIVo . (2.13) 

The brane is located dX y = y^ such that 6(?/o) — 1? that is, 

sinh(yo/^) = . 
In this geometry, the equation of motion for the scalar field is 

Yidy[b^dy<p) - ^ (0 + 3i7(/) - ^ aV) = . (2.14) 
This equations is separable and one can solve it by looking for a solution of the form (f) = f{y)(p{t, a:*), with 

if + SHif ^A<^ + mV = 0. (2.15) 

The separation constant corresponds to the square of the KK mass, as measured by an observer on the brane. 

Since there is no coupling between the brane and the bulk scalar field, the boundary condition for the scalar field 
at the brane location is simply dy(j) = 0, and therefore dyf = 0. The equation along the j/-direction implies that the 
mass spectrum is characterized by a mass gap 3H/2 p^. The corresponding eigenfunctions / can be written in terms 
of the associated Legendre functions. 

Let us now focus on a single KK mode, which is spatially homogeneous and sufficiently massive: m 3> H. One 
finds from I2.15|l 

fi*) = -4h cos(mi) . (2.16) 

If we take a time average over a time scale much longer than the period of oscillation m~^, we can ignore the oscillatory 
behavior and use 

(sm^imt)) = (cos^(TOi)) = i , etc. (2.17) 



From Eq. H2.8() , we thus find 



a3 



^(''^ = ^1/™!'"^'^'^], (2.18) 

where is the value of f{y) on the brane for the eigenvalue m^. From Eqs. (|2.1U|) and H2.16|l . and from the fact that 
dy(f> — —m?if) on the brane, we can evaluate E^^, as 

-E^,^-^^Kl\fm\^m^^5]^ (2.19) 
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where we have neglected the terms that depend on the initial data, which behave as and thus become negligible 
at late times. 

The above results show that the Weyl term Efj_,y contributes negatively to the effective energy density and pressure 
on the brane for a massive mode. Moreover, if one computes the total effective contribution of the bulk, i.e., the sum 
of T^^'^Q,^ and of the Weyl term Ea/s, one finds for the effective energy density and pressure on the brane 

2 / 2rr{b)t Tpt \ 1 2|j- |2 2 1 

KiP(cS) = - [k^T'^ > t-E tj = --Ksl/ml m — , 

4P(eS) = I {kIt'^''>\ - = . (2.20) 

This represents the backreaction effects of the bulk scalar field, which are of order 0(0^). Whereas the effective 
pressure due to the KK mode vanishes, because the bulk component and the Weyl component exactly cancel each 
other, the effective energy, remarkably, is negative. 



C. KK mode for a low energy cosmological brane 

We now calculate the effective energy density and pressure of a KK mode for a low energy cosmological brane. The 
bulk geometry around a brane, with a flat FLRW geometry and located at y = is given by the metric 0, 0| 

ds^ = -N'^it, y)dt^ + Q'^{t, y)a^{t)dx^ + dy'^ , (2.21) 

where 

Q(i,2/) =cosh(2//^)-77sinh(|y|/£) 

N{t, y) = cosh(y/^) sinh(|2/|/^), (2.22) 



H 



with 



y/mP + 1. (2.23) 



We have assumed that there is no dark radiation, i.e., that the bulk geometry is strictly AdS and not Schwarzschild- 
AdS. In general, this metric is non-separable. However, in the low energy limit characterized by HI <C 1 and Ht^ ^ 1, 
we have ~ 1 and f]/ H <C 1 so that the metric can be approximated by 

which is now separable. If one considers the evolution of a massless, minimally coupled scalar field in the above 
background metric, one finds that the field equation is separable and thus admits a solution of the form (j){t^ y) — 
f{y)(p{t) with 

<f + 3H<f + m^ip = , (2.25) 

where the function f(y) is assumed to be Z2-symmetric. 

The solution for f{y) with the appropriate Neumann boundary condition on the brane, /'(O) = is given in terms 
of the Hankel functions. There is a zero mode corresponding to m = as well as a continuum of KK modes with 
771 > 0. For a massive KK mode m ^ H, the four-dimensional part evolves as 

f = -^j7^ cos{mt) , (2.26) 

Similarly to the de Sitter brane case, one can compute the projection of the bulk energy-momentum tensor on the 
brane and one finds for its components: 
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This gives 



I to 

and £*t = — Thus we obtain 



r^*")', - Kg^i;', = 0, (2.29) 
at late times. Therefore, the effective energy density and pressure for a KK mode becomes 

«lP(eff) = \ - E\) = . (2.30) 

This means that, also for a low energy cosmological brane, a massive KK mode behaves as cosmic dust with negative 
energy density. 

The analyses given above imply that the result is independent of the existence of a mass gap and the essential 
factor is the background expansion of the brane. A KK mode can be approximately defined only for a cosmological 
brane which slightly deviates from the dS geometry and for a low energy brane, thus we expect that our result can 
be applied at least for these cases. However, for intermediate energy scales a KK mode is not well-defined in general 
and it is not clear how our result might be applied. 

Finally, we note that the bulk energy density of a KK mode on the brane remains positive as 

i4p{huXV) ■■= -i4T\ = \i4\fm?rn^^ > , (2.31) 

for both de Sitter and low energy branes (with the understanding that the time average over scales greater than 
is taken). It shows that there is no singular effect in the bulk in contrast to the peculiar behavior on the brane. 

III. EFFECTIVE THEORY IN THE BULK AND ON THE BRANE INCLUDING THE 

GRAVITATIONAL BACKREACTION 

After having studied the backreaction of the KK modes of a bulk scalar field, we now turn to the main subject of 
this paper, which is to study the backreaction of the gravitational perturbations of the metric itself on the cosmology 
of the brane. In this section, we adopt a more general perspective by considering a [d — l)-branc embedded in a 
(d + l)-dimensional bulk spacetime, although we remain primarily interested by the case d = A. This allows us to 
investigate the dependence on the number of dimensions of various quantities introduced in this section. 

A. Effective theory in the bulk 

We now consider only pure gravity in the bulk. The action of the system is given by 

S[9] = J d^+^x^i^^^+^^R - 2Ad+i) - j d^x^qa, (3.1) 

where A^^+i is the bulk cosmological constant and a is the brane tension. We mainly consider a dS brane background 
in this section and assume that its tension is larger than that of the corresponding RS value 2(d— 1) / {k^_^-^1) , where 
I = {-d{d - l)/(2Ad+i))V2 is the bulk AdS curvature radius. 

(0) 

We start from an unperturbed metric 9 , which is a solution of Einstein's equations and thus satisfies 

^[5](o)=0, (3.2) 
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where and in what follows the notation, Qfa + , means that a functional Q[a + g] of g is evaluated for a function 

/ 

/, i-e., 



Q[a + g] ^^Q[a + f]. 



(3.3) 



(1) 

We then consider (small) linear perturbations of this metric, which we write e g and such that its average vanishes 
i.e., 



(.9) 



0. 



(3.4) 



. (1) . 

Here we should specify our definition of averaging. We assume that the perturbation g has a typical wavelength A 

(0) 

which is much smaller than the characteristic curvature radius L of the background g , X L. Then we take the 
average over a length scale much larger than A but much smaller than L. In our case, we can take this average in the 
spacetime dimensions parallel to the brane. However, the situation is dramatically different in the direction of the 
extra spatial dimension because the brane is infinitesimally thin, which implies that the curvature radius along the 
extra dimension is infinitely small. Therefore one cannot take an average in that direction at or around the brane. 
Thus our averaging will include only the average over the 1 + (d— 1) spacetime dimensions. (For spatially homogeneous 
perturbations, we take only the time average.) 

What we are interested in is the correction to the original metric due to the backreaction of the metric perturbations. 
The total metric we consider can thus be written as 



(0) 



(1) 



(2) 



5tot ^9 +e g 9 



(3.5) 



(2) 

where the quantity 9 represents the backreaction due to the metric perturbations, so that the effective background 
(homogeneous) metric, after averaging, is given by 



g 



(0) 

= g 



2 

-e^ g 



(3.6) 



For convenience, the parameter e is introduced as an expansion parameter, which is to be set to unity at the end of 
the calculation. 

(1) 

If we expand the action with respect to g , we have 



r (1)1 


= S 










g + e g 


9 


+ 7fW 







(3.7) 



(1) 

Hence the variation of the above expression with respect to 9 yields 

SS 



SS, 
og 



(1) 

e g 



SS 
og 



6g 



(3.8) 



where we have used Eq. H3.2|l in the final equality. This implies that, up to 0{e), the equation of motion for the 
perturbation .9 is given by 



SS AO) , 



(1) 



= 0. 



(3.9) 



On the other hand, the variation of the action with respect to gtot gives 

gtot Sg 



~g 

5S_ 
Tg 



(1) 

a 5g^ ^ ' 
a 5g^ ^ ' 



( 



(i; 

e .9 



(1)\2 
(1)\2 



(3.10) 
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(0) 

where, to get the last expression, the argument of the coefficient of the third term, g, has been replaced by 9 , which 
is justified within the accuracy of 0{e^). If one averages the above expression, the second term on the right-hand side 
vanishes and we obtain the equation that determines the backreaction-corrected background metric g, in the form 



1 2/(1) S^S 



(1) 



(3.11) 



Substituting the explicit form for the braneworld action, we find that Eq. (|3.11|) yields 

(d+l)Qa^ + Ad+lS'^b = l^\+lT°'b + ^(braiic)°6 + <5i (branc) °fc , (3-12) 

where ('^+i)G is the background bulk Einstein tensor including the backreaction effects, i.e., for the metric g. And 
the stress-energy tensor due to the backreaction in the bulk is given by 

«^+ir%^-(('^+i)G\), (3.13) 

(2) 

where Q°■^^ is the bulk Einstein tensor at quadratic order. Here it may be worth noting that averaging is necessary 
for this effective stress-energy tensor to be physically meaningful, since there exists no locally covariant gravitational 
energy-momentum tensor due to the equivalence principle. The tensor t(brano)°6 corresponds to the brane energy- 
momentum tensor in the background configuration defined by the metric g and thus comes from the variation of the 
brane action in the left-hand side of H3.11|l . Finally, <5i(brane)"bj which comes from the brane-dependent part in the 
right-hand side of H3.1HI . denotes the backreaction due to the brane fluctuations and will be discussed in Section III. C. 
The existence of this term is the most important difference when compared to the case of the scalar field, in which 
case the backreaction originates purely from the bulk. 

Hereafter, we write ('*+i)G°b as ('*+i)G°b for simplicity. For the moment, we concentrate on the effective theory in 
the bulk, 

(^+i)G% + Kd+i5\ = «d+ir% . (3.14) 

Our first task is to evaluate the effective bulk energy-momentum tensor T°f,, which is quadratic in the metric pertur- 
bations. Then we will take the limit to the brane. 

We now identify the background metric g'"' with the separable metric of AdSd+i bulk-dS brane spacetime and g'l) 
as the linear perturbation of this system. Namely, 

ds" = dy^ + b'{y){l^.. + h^,)dx^'dx\ /i"o = /i/ 1/3 = 0, (3.15) 

where h{y) is the warp factor defined in Eq. (|2.12() and 7^^ is the metric of a d-dimensional dS spacetime which is an 
extension of Eq. H2.13|l . Note that we have adopted the so-called RS gauge for the perturbations 0, . The equation 
of motion for the perturbations in the bulk reads 



Kp^O. (3.16) 



This equation is separable and one considers solutions of the form hap — f{y)Lpap{x^)^ where f{y) is the generalization 
of the solution of Eq. (|2.15|l to the case of a d-dimensional brane with boundary condition dyf{y) —Qaty — yo because 
dyha/3 = on the brane. Similarly to the scalar case, the separation constant m represents the effective mass of a KK 
graviton mode and satisfies m > {d — l)H/2. The d-dimensional part ip"'^ satisfies 



¥>"^=mV"/3. (3.17) 



^ H^. Furthermore, for simplicity, we focus on perturbations of the tensor-type 
with respect to the spatial (d — l)-geometry, namely on those with h*t = h*i = h^t = 0. Taking the slicing of the de 
Sitter space with the flat spatial (d — l)-geometry, they will have the form, 

^^ = ^(^T)72 ' (3-18) 

where fm is the amplitude of the KK mode and Q^j is the polarization tensor on the flat (d— l)-space. The amplitude 
fm can be determined, for instance, by the normalization condition if one considers a quantized perturbation theory. 
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As mentioned earlier, in order to obtain the stress-energy tensor that embodies the backreaction due to the metric 
perturbations, one needs to "average" the Einstein tensor at quadratic order, according to Eq. H3.13|l . The components 
of the bulk curvature tensors, up to quadratic order in the perturbations are listed in Appendix A. As explained after 
Eq. I|3.4|) . we take the spacetime average in the 1 + (d — 1) dimensions parallel to the brane, but not along the extra 
dimension. In particular, because of the cosmological symmetry, we can take the average in the {d — I) dimensions 
over the complete space. The derivatives along the extra dimension are replaced by using the field equation (|3.16(l 
and the boundary conditions on the brane. Our procedure is detailed in Appendix B. 

Using Eq. ljA7|l of Appendix A and the computational rules detailed in Appendix B, we obtain in the limit y +0 
the expressions 

(''+'^ G = -^(/i-^D^v) - ^S^^h^^D^hp^) - l^/,p-l"/,^^l^^ . (3.19) 

A priori, the effective energy-momentum tensor includes an anisotropic stress, to which each mode will contribute 
with a factor 0{m?). However, if the perturbations are described by a random field which is statistically homogeneous 
and isotropic, the average over all modes of the anisotropic part must cancel. What remains is thus to justify the 
randomness of the perturbations. In this respect, the quantum fluctuations are indeed expected to have this property. 

(2) 

Also, (('^+1) u) vanishes on the brane by using the boundary conditions dyhap = on it. 



B. Backreaction on the brane 



Let us now discuss the effect of the backreaction onto the brane. The projected gravitational equation on the brane 
reads 



where 



A, 



d-2 



eff 



-A 



d-2 



d+l 



(3.20) 



(3.21) 



is the effective cosmological constant on the brane, and 



P = 



d-2 



d-l 

d-2 
2{d-l) 



(2) 



G 

(rf-2)(d-3),„ 



(2) 



_ ±/{d+l) Y, a 
d\ ^ ° 



M{d-\) 



pa 



d-2 



pa\f3 



(3.22) 



is the projection of the effective energy-momentum tensor of the bulk gravitons. The tensor t"^, corresponding to 



5t, 



(bra 



of the previous subsection, describes the brane perturbation induced by the bulk perturbation. We will 



show in the next subsection that, for our purposes, this term can be neglected. We now concentrate on the effect of 
the effective energy-momentum of the bulk gravitons projected on the brane, i.e., the terms T'*"^"^ and i?"^. 

Let us first consider T^^'^'^/j. Because of the assumed symmetries, i.e., the spatial homogeneity and isotropy, this gives 
in the brane an effective perfect fluid with some energy density and pressure. Decomposing the metric perturbations 
into KK modes, one finds that the contribution of a sufficiently massive mode to the energy density and pressure is 
given by 



{d + 3){d-2) 1 / 
16d(d-l) a-i-i"" ^ 
{d^ + 3){d-2 



kl / 1 



lM{d-l) 



(3.23) 



We must also take into account the projection of the Weyl tensor on the brane, E°- p. Although this term is not 
included in the "bulk energy- monemtum" tensor because it is a part of the bulk Weyl tensor, it contributes nevertheless 
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to the projected gravitational equations as an "energy-momentum" tensor. Although its direct evaluation is rather 
delicate, this term can be computed by resorting once more to the cosmological symmetry. From Eq. (I3.20|) . the 
contracted Bianchi identities D°'^'^'^Ga^ = 0, together with the conservation of Tq,^, give 

D^'E^, = kI+,D^T(^J. (3.24) 

Because of the cosmological symmetry, the only non-trivial component of the above equation is the time component, 
which reads 

d,E\ + d^E\ = «2 UTi^)t^ ^ _ 1)^tO^)*, - ^rC^)^) , (3.25) 
a \ a a J 

where, on the left-hand side, we have used the property that E^i, is traceless and thus E'^i = —E^f The integration 
then yields 



^ t dt'a'ldtT^'^^'t + (d- l)^r('^)*, - ^T^^)\) . (3.26) 



As before, we neglect the contribution from the initial condition, which is valid at late times. 

Substituting a KK graviton mode given by Eq. H3.18() into the integrand on the right-hand side of Eq. H3.26|l . and 
taking the time average, one finds 

This gives, at late times. 



Because of the traceless nature of this tensor, we then obtain E^j — —{l/{d~ l))E*tSj. 
The total contribution of the two tensors is therefore 

^2^^y(b).^ _Et^^d-^J_ m^\f^\^(Q>^^Ql,) . (3.29) 

for the temporal part and 

'^l+lT^''^'^ - E\ = 0, (3.30) 

for the spatial part. This means that the contributions of a KK mode to the total effective energy density and pressure 
are respectively given by 

2 



For instance, for d = A, we obtain 



KdP(off)=0. (3.31) 

^lP{.n) = . (3.32) 

The effective isotropic pressure vanishes and the effective energy density is negative. This is the same as in the case 
of the scalar field discussed in the previous section. 

We note that the bulk energy density of a KK mode on the brane remains positive 

'«^+iP(buik) := = ^^,^''\U^{Q'''Qm) > 0, (3.33) 



a" 



as in the scalar case, Eq. (|2.31|l . It shows again that there is no singular effect in the bulk. The negativity of the 
effective energy density on the brane originates from the projected Weyl tensor E^n. 
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C. Brane intrinsic contributions 



We now consider the branc intrinsic contributions. In order to discuss the gravitational perturbations in the brane 
world, it is not sufficient to consider the contribution from the bulk. The brane perturbations must be taken into 
account as well. We take an approach in which we derive the second order boundary action and regard it as the action 
for an effective matter on the brane. 

In this paper, the brane is treated as a thin wall. In the thin wall approximation, the second order action on the 
boundary has been derived in the Appendix of 20ij. When there is no ordinary matter on the brane and thus no 
brane bending mode, the second order boundary action is given by 20] 



(3.34) 



where /i^i/ — b^h^^, kp^r — dyhpcr/2 and AQ = Q^^-* — Q*- For an AdS-bulk configuration and with the assumption 
of Z2 symmetry about the brane, this reduces to 



— -a d'^xJ^hp^hP" . (3.35) 
4(rf-l) JoM 

The second order action can be regarded as an action for some effective matter induced on the brane 



3 

,,r^' 2(d- 1) 



d'^x./^Cm, An TTTT^^f^ V/i"" ■ (3-36) 



Its variation with respect to the background metric q^i, yields the induced matter energy-momentum tensor on the 
brane 



" d^'^i'^''"''''^^ ~ i^a/j/j^'V) • (3-37) 



Ta[3 



Note that, strictly speaking, qp^, does not include the backreaction. However, as discussed in Section III. A, for 
perturbations with small amplitude, the linear perturbation equations are identical to those for the background 
metric in which the backreaction is taken into account. Thus we can add this term as a part of the (effective) matter 
contribution in the effective equation on the brane. 

We can readily calculate the effective energy density and pressure of this contribution. One finds 

_ 3(rf — 2) ^4 2\p |2 1 /r)ke^ 

16{d- 1)2 f^+i I-'™! Tfl^^X*^ ^ke/^ 



K<iP(brano) - 16(d - 1)2 '^d+l'^ l-^™' Q*k 

3(d-2)(-d + 5) 4 2,. ,2 1 /r^kir.* \ 



KdPihianc)^ lg(^_ 1)3 '^d+l<^ ^d^X^ Q ki / ■ (3.38) 

This gives the equation of state 

W^Cbranc) ^'JZ^ ' (3.39) 

For d = 4 the boundary contribution thus behaves as radiation, for d = 5 as dust and for d > 6 as matter with 
negative pressure. However, its contribution to ^''-'G"/3 is of order 

where we use |2r)( 

2(d-l)/ ,,,,^2\i/2 



a = 



Hi)) . (3.41) 



For the cases Hi <t: I and Hi > 1, the right-hand side of Eq. IpHHIl is 0{i-^h'^) and 0{H^h'^), respectively. Thus 
as long as we consider sufficiently massive KK modes, with m ^ max{^~^, H}, the brane perturbations can be safely 
neglected and only the projected bulk contributions are relevant for the effective theory on the brane. 
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IV. NEGATIVE ENERGY DENSITY FROM THE BULK POINT OF VIEW 



Intuitively, the "negative energy density" of a KK mode is rather puzzhng. However, we can understand its cause 
by regarding the KK modes as a part of the dark component, such as the dark radiation. The energy density of the 
dark component evolves as 

P(D) + 4.Hp(^D) = -2(1 + ^)Tafcu"n'' - 2{He)Tabn''n'' , (4.1) 

where and n° are the tangent and normal vectors to the brane, respectively. Since there is no matter on the brane 
and no brane-bulk energy exchange, the first term on the right-hand side of Eq. H4.1|l vanishes and only the second 
one, related to the pressure transverse to the brane, survives. In terms of the energy conservation law in the bulk, this 
has the simple interpretation that the work done by the pressure on the brane to move it outward in the direction of 
the AdS infinity reduces the energy in the bulk. As a result, the dark energy density decreases, since the dark energy 
density on the brane is proportional to the total mass (energy) in the bulk |0| . For a massless scalar field, 

Ta^nV = 102-^0. (4.2) 



Thus, the dark component decays faster than ordinary radiation. For the KK modes, after time averaging, we have 

1 

4^ 



T„fc7iV = -i^|/,„|W. (4.3) 



The formal solution of Eq. (|4.1|l is 



'to 

where C denotes the initial mass in the bulk. For the KK modes, 



PiD)=-^ f dtaHHe)Tabn'^n'' + ^, (4.4) 



f dta\H£)Tabn-n'' = Uu^m-" f dt d ^ j (a{t) - a{to))\frn\^m\ (4.5) 

JtQ ^ 'J io 



Hence 

PiD)=-^\Wm' + ^, (4-6) 



where we have redefined the mass parameter C by absorbing into it the initial data dependent term of the integral 
(ESI). 

Anyway, in the case of a dS brane (or a cosmological brane which slightly deviates from the dS geometry), the 
effective cosmological constant dominates the cosmological evolution and the KK effect does not have a significant 
impact on the brane. For a low energy brane, especially for a radiation-dominated brane, naively one might worry 
that this result would imply the appearance of a negative energy density within a finite time. However from Eq. (|4.1|l . 
the bulk pressure term is proportional to H. Hence ii H < a,t H = 0, the energy density will remain positive at the 
expense of rendering the universe to recoUapse. 

For simplicity, we consider the case where the cosmological evolution of the brane is determined solely by the dark 
component. Note that this discussion can be generalized when one considers ordinary dust or radiation in addition 
to the dark component. The Hubble parameter on the brane is obtained from 

= ^p,^. = + — 4 = -^rabn'^n" + ^^. (4.7) 

Taking the time derivative of this equation, we obtain 



- f - If ? = If ^- <"> 

Therefore, at H = 0, we have 

1 



= -^r.bnV <0, (4.9) 
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and the universe begins to collapse. Thus, the backreaction of the KK modes leads to a collapsing universe. 

The situation is the same for the case of KK gravitons as long as the brane fluctuations are negligible, because we 
have 

nlTabu'-n'' = -( G \) = im^^/i^^v) = ^[/..^^^^(g^^QI,,) > 0. (4.10) 

Thus, provided that the brane fluctuations can be neglected, the brane universe will start to collapse within a finite 
time. For more realistic situations in cosmology, our result suggests that for a low energy brane the brane universe 
will eventually collapse unless the contribution of the true (normal) dust matter is larger than that of KK modes. 

V. SUMMARY AND DISCUSSION 

We have investigated the effect of a Kaluza-Klein (KK) mode on brane cosmology, focusing on the second Randall- 
Sundrum (RSII) type model, i.e., for a Z2-symmetric brane embedded in the anti-de Sitter bulk. 

The KK gravitons, which are just the metric perturbations in the bulk, are produced during a de Sitter (dS) brane 
inflation phase via vacuum fluctuations. From the four-dimensional point of view they are effectively equivalent to 
massive gravitons with masses m > 3-ff/2, where H represents the dS expansion rate of the brane. The theory of 
linear perturbations reveals that the squared amplitude of a KK mode decays as a^^ and its contribution rapidly 
becomes negligible during the brane inflation. However, after brane inflation in the radiation-dominated Friedmann- 
Lemaitre- Robertson- Walker (FLRW) era the background radiation energy density decays as a^^, which implies that 
the contribution of KK gravitons may have a significant impact on the brane cosmology at late times. 

Before discussing the case of the KK gravitons, we have considered the case of a massless, minimally coupled 
bulk scalar field as an exercise. This is because a massless, minimally coupled scalar field has many properties in 
common with the gravitational perturbation, but it is much simpler to deal with the former than with the latter. We 
have considered two limiting cases of the background spacetime in which the equations of motion become separable, 
namely, the case of a de Sitter brane and the case of a low energy cosmological brane. In both cases, we found that a 
sufficiently massive KK mode, with mass much greater than the Hubble expansion rate, H, indeed behaves like 
dust but, rather surprisingly, its energy density is negative. 

Then we have turned to the case of the bulk gravitational perturbations. We have first derived the effective energy- 
monemtum tensor of a KK mode in the bulk, and investigated the effect of the KK mode on the brane cosmology by 
projecting the effective energy-monemtum tensor on the brane. We have found exactly the same behavior as that of 
the scalar, i.e., a massive KK mode behaves as a negative energy density dust. 

The negative energy density of a KK mode may sound rather puzzling. But, from the bulk point of view, we have 
shown that this result can be regarded as a natural consequence of the energy conservation law in the bulk. Here 
the essence is to recall that the so-called dark radiation term, which behaves like radiation on the brane, describes 
the total mass in the bulk. Then, a very massive KK mode corresponds to a particle with a high momentum in the 
direction of the extra dimension, which exerts a pressure on the brane and pushes it outward in the direction of the 
AdS infinity. As a result, the energy in the bulk decreases, leading to the decrease of the dark energy term. Thus, a 
massive KK mode gives a negative contribution to the dark radiation term. This is why a KK mode behaves like a 
negative energy density dust. 

Note that the negative energy of a KK mode emerges only from the effective four-dimensional point of view on the 
brane. The bulk energy density for a KK mode still remains positive and thus there is no singular effect in the bulk. 

We have studied the two cases (de Sitter and low energy branes) in which the bulk equations are separable, hence 
the KK modes are well defined. However, for a general cosmological brane, one cannot define a KK mode since its 
very definition depends on the separability of the equations in the bulk. Nevertheless, considering the discussion from 
the bulk point of view given in the previous paragraph, it seems reasonable to expect that this backreaction effect 
of the bulk metric perturbations persists for a general cosmological brane. Thus we conclude that the effect of very 
massive KK modes is to reduce the energy density on the brane, hence the expansion rate, and for a low energy 
brane the universe will recoUapse unless the contribution of a normal (true) dust matter is larger than that of the KK 
modes. 

To quantify this effect in realistic cosmological models, there are some additional issues that remain to be resolved. 
In this paper, we have considered only a single KK mode and calculated the effective energy density and pressure. In 
reality, one should integrate over all the KK modes that contribute to the cosmology of the brane. This requires first 
the knowled ge o f the whole spectrum of the KK modes, which will presumably be determined by vacuum fluctuations 
in the bulk (2lL 122 . l2^. However, knowing the whole spectrum may not be enough, because a naive integral of the 
KK spectrum is expected to diverge. One would then need an appropriate regularization scheme. In connection with 
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this, it may be important to take into account the thickness of a brane. either classically as in the case of a classical 
domain wall or quantum mechanically by considering the wavefunction of a brane. 
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APPENDIX A: SECOND ORDER CURVATURE TENSORS IN THE BULK 

Here we spell out the components of the curvature tensors up to quadratic order in the bulk metric perturbation. 
We consider the {d + l)-dimensional perturbed metric in the form, 

ds^ = dy^ + b\y) (7^, + /i^,) dx^'dx'' , (Al) 

where 7^1^ is the metric of the background d-dimensional spacetime section. In the text, we identify 7^1/ with the 
metric of a de Sitter spacetime. We impose the following gauge conditions on the perturbation: 

= 1/3 = 0, (A2) 

where the vertical bar ( | ) denotes the covariant derivative associated with the d-dimensional metric 7^1/ , and the 
tensor indices of hfj,^, are raised or lowered by the metric 7^,^ (not by the 5-dimensional metric). 
The non-trivial components of the connection are given by 

h' 1 1 



where the prime ( ') denotes the t/-derivative. 

The non-trivial components of the Riemann tensor are given by 



{d+i)py =-lh'^(h , -h ,\-—h' (hP I +hP , -h 



b'- / 

_l h' hP , + hP , - h IP 



2 
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-hpJPhP^\l3 - V'"^"/?!" + hpJPhco^P - h^p^phP^i^ - /i%|/3/i'',.|a + h^piph^^^^P 

-hPp\phP^\, - h%^phP,\^ + hPp\phJp + ftp/jl^ft^i, + hp0^PhP,^„ - hpp\^hjp) . (A4) 

The mixed components of the Ricci tensor are given by 

(c^+Dij.^ = -d^ + \hP^%^ + ^-hP^h'p^ + lhP^h% , 

('^+^)Ry, = -^hP'^K.l. + -hP^l.K'^ + \hP^K^\. 



1 lb' 1 lb' 

-^h^P'h'pf, + --S-0hP^h'p, + -h'^Ph"p0 + -d-i 



1l3p\a 



-^h-'iK'^l^^ + h'^lp^ - °d hpff) + l('^)R^0h-Php'y . (A5) 
The Ricci scalar is given by 

= -2d'^ + IC^)/? - d{d - 1) - ^^^^R,.hP^ 

+\hp^'h'p, + {d + l)^-hp''h'p, + hp^h';,, + ^hp^a^hp, 

+ ^ (i^"'^'" ViM - ^^'"^''^ Vk) - ^h^PW.hP^ + ^'■''^Rp'^ho.h'^'^ . (A6) 
Using these results, the components of the Einstein tensor are given by 



-IhP-^'h'p, - ^{d- l)^-hP^h'p^ - ^hP^a.hp, 



(^^+1)0"^ = {d- 1)^<5"^ + 1 WG"^ + i(d - l)(d - 2) (^) \ 



-^(^^Rp^h'^P + ±^id)Rp^hP'^5'-p 

+ lh"P'h'pp + \h'^''h"pf, + \d'^h'^%fi - ^hP^(hpp\-\„ + h%\p„ - h%\p„ - 
+ 4^ - 2/i"''l'^V|, + 2h'^P\'^hpp\,) + l(<^)Rp0h'^-h,P 
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-{- 



"'pa\ii 2 PPW 



4" "p' 



^2 pkM"- 



('i)R hP /iMCT 

^2 ^^pa'^ ^"^ 



(A7) 



APPENDIX B: COMPUTATIONAL RULES FOR AVERAGING 

In this appendix, we describe the computational rules for averaging the components of the second order part of 
the curvature tensors listed in Appendix A. As we have noted in the main text, the notation (A) includes both 
the averaging along the ordinary spatial dimensions which are assumed to be homogeneous and isotropic, and the 
small-scale time averaging as defined in H2.17II . In both cases, the computational rules are similar (See e.g., Ref. l23|). 
However, we do not apply the same rules for terms with derivatives in the bulk direction, because we are dealing with 
a braneworld and the averaging along the bulk direction is ill-defined. 

First, we note that we are interested in massive KK modes. So, we can neglect terms coupled to the background 
curvature tensor as 



(Bl) 



which are of order 0{h? / L'^), where L is the d-dimensional characteristic background curvature radius, in comparison 
with terms as 



lhf"^\ h I \ Ihf^h I \ 



(B2) 

which are of order 0{m?h'^). For instance for a cosmological brane with expansion rate i?, we have L ~ 0{1/H). 
Thus H implies m 3> L^^, and we can safely neglect corrections of the form HIjl(l . 

As a consequence, when taking the average, we are allowed to freely interchange the order of the covariant derivatives. 
For example. 



Ifii/hpcr^ — \vphpiy^ 



(B3) 



where corrections of order 0(}t? jl?) are neglected. From now on, as in the main text we will use instead of "~" 
by neglecting the corrections. 

Another computation rule is that total derivative terms can be neglected. For example. 



h'"'\ h 



hP^h 



pa\v 



hP^h 



'p(T\iyfj. 



hP^h 



(B4) 



This is because the total derivative term can be cast into the surface integral which is smaller in magnitude than the 
volume term by a factor mi? (3> 1), where R is the length scale of the averaging volume which is taken to satisfy 
R > m~^. 

As mentioned above, we do not apply the same rules for terms with derivatives with respect to the bulk coordinate 
y. However, when one considers projections onto the brane, some simplifications occur. On the brane, we have the 
boundary condition branc = 0, which enables us to neglect all the first derivative terms, e.g.. 



V/3 



ipp 



= 0. 



In addition, using the bulk equation of motion (|3.16() . we have 

'hP'' Dd hp„ 



branc 



branc 



(B5) 



(B6) 
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